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ABSTRACT 

The question of determining the spatial geometry of the Universe is of greater 
relevance than ever, as precision cosmology promises to verify inflationary predic- 
tions about the curvature of the Universe. We revisit the question of what can be 
learnt about the spatial geometry of the Universe from the perspective of a three-way 
Bayesian model comparison. By considering two classes of phenomenological priors 
for the curvature parameter we show that, given current data, the probability that 
the Universe is spatially infinite lies between 67% and 98%, depending on the choice 
of priors. For the strongest prior choice, we find odds of order 50:1 (200:1) in favour 
of a flat Universe when compared with a closed (open) model. We also report a ro- 
bust, prior-independent lower limit to the number of Hubble spheres in the Universe, 
Njj > 5 (at 99% confidence) . We forecast the accuracy with which future CMB and 
BAO observations will be able to constrain curvature, finding that a cosmic variance- 
limited CMB experiment together with an SKA-like BAO observation will constrain 
curvature independently of the equation of state of dark energy with a precision of 
about (7 ^ 4.5 X ID""*. We demonstrate that the risk of 'model confusion' (i.e., wrongly 
favouring a flat Universe in the presence of curvature) is much larger than might be 
assumed from parameter error forecasts for future probes. We argue that a ba detec- 
tion threshold guarantees a confusion- and ambiguity-free model selection. Together 
with inflationary arguments, this implies that the geometry of the Universe is not 
knowable if the value of the curvature parameter is below jri^l ^ 10^^. This bound 
is one order of magnitude larger than what one would naively expect from the size of 
curvature perturbations, 10^^. 

Key words: cosmology: theory; cosmology: cosmological parameters; methods: sta- 
tistical. 



^ 1 INTRODUCTION 

Constraints on the total energy density of the Universe, 
Sltot have improved spectacularly in the last two decades. 
Before the onset of precision cosmology, the total matter 
energy content of the Universe was known only with order- 
of-magnitude precision. The determination of the angular 
scale of the first acoustic peak in the cosmic microwave 
background (CMB) was a major milestone towards deter- 
mining the spatial curvature. The location of the first peak, 
£ ~ 220, together with estimates of the Hubble constant, 
implies that the Universe is close to flat. While 10 years 



ago this s tatement could be made with an accuracy of or- 
der 10% (jde Bernardis et aLllioOoh . more refined measure- 
ments of the CMB power spectrum by WMAP and other 
experiments have reduced the stat istical uncertainty to s ub- 
percent precision in recent years (jKomatsu et al ]|200i). In 
turn, this has allowed us to tighten constraints around a fiat 
Universe with no spatial curvature, S7k = 1 — Qtot ~ 0. This 
spectacular increase by over a factor of 100 in accuracy in 
less than two decades refiects huge steps forward in detector 
technology, telescope design and computing power. 

As there are only three discrete possibilities for 
the underlying geometry in a Friedmann-Robertson- 
Walker Universqj (namely, fiat, close and open, see 
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^ Although the space of models could be extended to include 
non-trivial topologies, in this paper we shall keep v^ith the sim- 
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however iMersini-Houghton et al.l (|2008l ) for a landscape- 
motivated alternative with an oscillatory curvature term), 
the question of which one of these three models is the cor- 
rect description for our Universe is particularly well suited to 
be phrased in terms of model comparison. In his pioneering 
application o f the Baye sian model comparison framework 
to cosmology, Ijaffd (fl996. ') found that the determination of 
the Hubble parameter using the Cepheid variables method 
coupled with a lower limit to the age of Universe already 
allowed one to infer that a Universe with vanishing curva- 
ture and non-zero cosmological constant was the preferred 
model (albeit only with modest odds of about 7:1). More 
recently, the WMAP 3-year data implied that the odds in 
favour of a flat U niverse increased to between 29:1 (jTrottal 
l2007al ) and 48:1 (jKunz et al.ll2006l l when comparing a flat 
Universe with curved models (both open and closed). 

Unlike many dark energy models that are mostly phe- 
nomenological, models predicting curvature are rooted in 
fairly well understood physics, a feature which helps in set- 
ting physically motivated priors on the model parameters. 
For example, the possibility of a flat, £7^ ~ Universe hcis 
long been favoured by theoretical prejudice, as a flat or 
close-to-flat Universe is a generic prediction of the infla- 
tionary scenario, which appears to have been confirmed by 
observations to date. With the prospect of even more vigor- 
ous observational campaigns coming up in the next decade, 
it is timely to ask to which point the accuracy in n„ can 
and should be pushed before the question about the flatness 
of the Universe becomes irrelevant, uninteresting or unde- 
cidable. Determining curvature is also important in order 
to avoid mistaking a non-flat Universe for an indication o f 



an evolving dark ener g y density, se e e.g. Knox et ahl (|2006h ; 
IClarkson et al] (|2007l ): Ivirev et al.l (|2008l v 

In this paper we address the capability of future CMB 
and baryonic acoustic oscillations (BAO) observations to 
constraint curvature, both from the point of view of param- 
eter constraints and from the perspective of a three-way 
Bayesian model comparison. We are primarily interested in 
the accuracy that can be achieved using the acoustic scale as 
a standard ruler, although complementary observations (for 
example, supernova type la (SNIa) or weak lensing observa- 
tions) will help to break existing degeneracies bet ween cur- 
vatur e and the dark energy equation of state (Clarkson et al.l 
I2OO7I ). thereby improving the statistical power. A fundamen- 
tal limit to our ability to determine curvature is set by the 
order of magnitude of local fluctuations in the spatial cur- 
vature, Af2tot ~ 10~^. We investigate how this translates 
into terms of model selection and, crucially, model confu- 
sion, and show that the size of the fluctuations means that 
the question of curvature becomes statistically undecidable 
for < 10~*, i.e. about one order of magnitude above 
the naive expectation, and this is regardless of the amount 
of data gathered. 

This paper is organized as follows. In section [2] we in- 
troduce the data we use and our forecast procedure, while 
we briefly review relevant aspects of Bayesian model selec- 
tion in section O where our prior choices are discussed. We 



plest option, namely that the Universe's topology is trivial, as 
searc hes for non-trivial topologies have been unsuccessful to 
date llCormsh et al.ll2004) . 



present the evidence from current data in section while 
section [S] gives the results of our forecast for future probes 
and discusses model confusion. We give our conclusions in 
section [S] 



2 SETUP AND METHODOLOGY 
2.1 Measuring the acoustic scale 

The acoustic peaks in the CMB power spectrum measure the 
projected sound horizon at recombination. The comoving 
sound horizon at decoupling is given by 



C 
Ho 



Hiz) 



dz, 



(1) 



where Hq is the Hubble constant today, z^cc is the redshift of 
decoupling and Ca is the sound speed of the coupled photon- 
baryon fluid. 



(2) 



\/3(l + i?)' 



with R = 3pt/pj ^ [670/(1 + z)]{Q.bh'^ /0.022). Here, pi 
and pj are the time-dependent energy densities of baryons 
and photons, respectively, while Qb is the energy density 
parameter for baryons today. The function H{z) is given by 



H^{z) = (n™(l + zf + nr{l + z)^ + Q.^{1 + zf 



+ f2dc exp ( 3 



+ w{x) 



+ x 



dx 



), (3) 



where the dark energy time evolution is described by 
the present-day dark energy density in units of the crit- 
ical density, fide, and by its equation of state, w{z). 
The energy density parameter for radiation (photons and 
neutrinos, taken here to be massless) is fir = ys'll + 
(21/8)(4/ll)''/3]T,iMB/^^ ~ 4.13 X IQ-^/h^, while 



alHl 



(4) 



is the curvature parameter (ao is the scale factor today). 
The curvature constant k determines the geometry of spatial 
sections: k = for a flat Universe, k = -1-1 for a closed 
Universe and k = — 1 for an open Universe. 

The comoving distance to an object at redshift z is given 

by 



dx 



Hoaa Jo H{x) ' 



(5) 



Given knowledge of the comoving length A of an object at 
redshift 2, measurement of the angle subtended by it on the 
sky, 9, determines its angular diameter distance, Da{z) 



Da(z) = 



A(a/ao) _ aogf.(x) 
9 ^ 1+z ' 



(6) 



where SK.{y) is y, sm{y) or si nh(?/) for k = 0, -1-1, —1, respec- 
tively . A number of authors (iBond et aL 19971: Bowen et al.l 



l20_o j ; i Melchiorri fc Gri flithj I2OOII : iKosowskv et al.1 l20oi 
I Jimenez et al.,, 2004 1 have pointed out that the morphology 
of the acoustic peaks in the CMB power spectrum is largely 
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controlled by the baryon density flth^ and by two 'shift pa- 
rameters' 



la = 7rx(zdcc)/rs(2doc), 



(7) 
(8) 



In the context of the recent interest in dark energy, the use- 
fulness of employing both shift parameters (and their corre- 
lations) as a handy summar y of CMB constraint s has b een 
brought into sharp focus bv lWang fc Mukheried (|2007l ). In 
this work, we follow their method of employing 'distance 
priors' as constraints on (la, R, Zdcc) for a summary of the 
information given by the CMB on the expansion history of 
the Universe. 

At lower redshift, the acoustic signatur e has been re- 
cently detected in the distribution of galaxies (jPercival et al.l 
I2OO7I : lEisenstein et al.ll2005l : ICaztanaga et al.ll2008[ l. thereby 
providing further constraints on the recent expansion history 
of the Universe. Future large galaxy surveys are expected to 
considerably improve present-day accuracy, by simultane- 
ously determining the angular diameter distance and the 
Hubble function H{z), which can be obtained by measuring 
the acoustic scale in the radial direction if spectroscopic data 
are available. This is because the radial extent of a feature 
along the line of sight is related to the redshift range Az by 



cAz 



(9) 



hence a measurement of ry allows a direct reconstruction of 
H{z). 

2.2 Parameters and data sets 

In this paper, we consider cosmologies containing baryons, 
cold dark matter, dark energy and a possible curvature term. 
The radiation density is fixed to the appropriate value for 
3 families of massless neutrinos throughout. Dark energy is 
taken to be either in the form of a cosmological constant, 
lu = — 1, or is described in terms of an effective equation 
of state WeS 7^ —1, which is taken to be constant with red- 
shift. Of course more complex parameterizations are possi- 
ble, and in particular one could consider an evolving dark en- 
ergy equation of state which cha nges with redshift ( for co n- 
straints on such models, see e.g. IZunckel fc Trottal l|200iD ). 
A particularly popular phenomenological parameterization 
of a time-evolving dark energy is to describe the equation 
of state as w{z) = wo + j^Wa, with two free parameters 
{wo,Wa)- We comment below on the impact that adopting 
such a more general dark energy model would have on our 
results. 

We employ a Metropolis-Hastings Markov Chain Monte 
Carlo procedure to derive the posterior distribution for the 
parameters in our model. We take flat priors on the following 
quantities: flmh^ , O^/i^, la, R, WcH (whenever the latter is not 
fixed to —1). The prior bounds on the first 3 parameters 
are irrelevant, as the posterior is well constrained within the 
prior. For WcS we take a prior range —2 ^ WcS ^ —1/3, with 
the lower bound cutting off some of the posterior for some 
of our data combinations (see below). Finally, the choice of 
prior for fl^ is fundamental for the model comparison part, 
and we discuss it in detail in section 

When considering present-day data, we include the 



WMAP 5-year data jPunklev et al.l |20 08l ) via their con- 
straints on the shift parameters and the baryon density , fol- 
lowing the method employed in iKomatsu et all (|2008E . We 
also include the SDSS baryonic acoustic scale measurement 
as an additional datum at redshift z — 0.35 by adding a 
Gaussian distributed measurement of the quantity 



A 



X (^bao) 



CZha, 



H{Zha. 



CZhs. 



(10) 



where Zbao = 0.35. We employ the mean value A = 0.47 4 
with standard deviation a a = 0.017 (|Eisenstein et al.]|2005l ). 
We also add the Hubble Key Project determination of the 
Hubble constant today, as a Gaussian datum with mean 
Ho = 72 km/s/Mpc and standard deviation ghq = 8 
km/s/Mpc (jPreedman et al.ll200L ). SN type la data are in- 
cluded i n the form of the S upernovae Legacy Survey (SNLS) 
sample l|Astier et all 120061 ). 



2.3 Future data 

We now turn to describe our procedure for simulating con- 
straints from future CMB and BAG observations. 



CMB data: Planck and CVL experiment 

We consider two types of future CMB measurements, one 
from the Planck satellite (due for launch early in 2009), 
which will measure the temperature power spectrum with 
cosmic variance accuracy up to ^ ~ 2000, and will consid- 
erably improve current precision in the ET and EE power 
spectra. We also consider a hypothetical Cosmic Variance 
Limited (CVL) experiment, which would measure the TT, 
EE and ET spectra with cosmic variance precision up to 
£ — 2000. This is meant to represent the ultimate precision 
obtainable from measurements of the acoustic scale at re- 
combination (although of course extra information on the 
expansion history will be available, e.g. via the integrated 
Sachs- Wolfe effect or CMB lensing. As mentioned above, we 
are concerned with the accuracy achievable by 'geometric' 
means alone). 

We start by choosing a fiducial value of the cosmolog- 
ical parameters around which to generate simulated CMB 
data. We employ Qth'^ = 0.02268, ficdm/i^ = 0.1081, fi^ = 
0,uioff = —1, which are in good agreement with the current 
best-fit from WMAP and other CMB observations (the val- 
ues of the spectral tilt and perturbations normalization are 
irrelevant for our analysis as we only employ effective dis- 
tance measures to the last scattering surface from the CMB). 
The correspondin g CMB power spectra are computed using 
the CAMB code fiLewis et al. 2000). W e then a dd no ise ac- 
cording to the procedure described in iLewid (|2005l ). with 
noise levels appropriate for either Planck or the CVL ex- 
periment (which has no noise up to ^ = 2000). Fi nally, a 
modified version of CosmoMC ()Lewis fc Bridlell200i ) IS em- 
ployed to fit the resulting noisy power spectra and to re- 
cover the covariance matrix for the parameters (R,l a , Zdec) , 
following the method described in iKomatsu et al.l (|2008h , 
which shows that constraints on this set of paramete rs are es- 
sentially equivalent to constraints on (R. L .Qhh^). \Li et al.l 
(2008) have analyzed in detail the loss of information in- 
volved in going from the full CMB data analysis to the use 
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^dcc 
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5.96-10-^ 1.96 ■10-"' 


1.02 ■ 10- 
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2.15 ■ 10^5 


1.06 ■ 10" 
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6.90 ■ 10- 


2 




Cosmic Variance Limited (CVL) 






la R 


^dcc 




la 


8.12-10-"' 3.89 10-5 


1.35 • 10- 
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6.23 ■ 10-*^ 


2.52 • 10- 
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^dcc 




1.47- 10- 


2 



Table 1. Covariance matrices for the distance parameters 
{la, R, z^gc) for Planck (top) and the CVL experiment (bottom). 



of the constraints on the set {R, la, Zdec) and have found that 
the covariance matrix method repr esents accurately t he in - 
formation contained in the CMB. iMukhcricc et all (j2008l ) 
investigated the application of this formalism to Planck pri- 
ors, and found a significant correlation between the shift 
parameters and the spectral tilt, ns- In this work, we do 
not include the tilt in the description of Planck data, on the 
basis that we never use Planck data alone to derive our con- 
straints on the curvature parameter. Thus, the degeneracy 
between ns and the shift parameters can be assumed to be 
effectively broken when including non-CMB observations, 
in particular data on the matter power spectrum which, by 
extending very considerably the lever arm of the CMB, are 
expected to be able to reduce the uncertainty on ns to a 
level which does no longer impact on the accuracy of the 
shift parameters. 

The fiducial values for our reference choice of param- 
eters are {R,la,Zdcc) = (302.06,1.709,1090.46). The corre- 
sponding covariance matrices for Planck and the CVL ex- 
periment are given in Table [T] 

In obtaining the covariance matrix for Planck and the 
CVL experiment, the curvature parameter has been allowed 
to vary (with a fiat prior over a suitably large range so that 
the posterior is much narrower than the prior), in order to 
obtain errors that correctly account for degeneracies in n^. 
On the other hand, the equation of state parameter has been 
fixed at w = — 1 when computing the covariance matrix. 
This is expected to be irrelevant as the whole point of using 
CMB 'distance priors' of this sort is precisely that they are 
largely independent on the assumed dark energy model (at 
least as long as the contribution of dark energy in the early 
Universe is negligible, which is the case here since we never 
fit evolving dark energy models). 

This covariance matrix is then used as the CMB high- 
redshift constraint. Notice that although the simulated data 
are obtained around a fiat fiducial model, we can safely use 
the resulting covariance matrix to represent CMB distance 
priors even when the fiducial model is slightly changed to 
Qn (as long as the change is not too large as to radically 
modify degeneracy directions), as we do below, where we 
employ fiducial models with ^ 5 x 10""^. 



BAO data: WFMOS and SKA-like experiment 

Regarding future BAO measurements, we adopt two bench- 
mark experiments. One is the Wide-Field Multi-Object 
Spectrograph (WFMOS), a proposed instrument for the 8- 



m Subaru telescope which will employ a fiber-fed spectro- 
graph to carry out a low [z ~ 1) and a deep (z ~ 3) survey 
to determine the acoustic oscillati on scale both in the trans- 
verse and in the radial direction (jBassett et ahllioosl ). WF- 
MOS could be operating around 2015. We also consider a 
more futuristic type of measurement, of the kind that could 
be delivered by the Square Kilometer Array (SKA) radiote- 
lescope around 2020 by performing a full-sky survey of HI 
emission. 

In modeling the accurac y of these o bser vations, we 
closely follow the treatment of iBlake et al. 

I (2006), to which 

we refer for full details. In summary, we employ the following 
fitting formula for the fractional accuracy of the determina- 
tion of the transverse and radial acoustic scale: 



X = X(i 



1 + 



off Djzof 
n blD{zf 



(11) 



with f(z) — (zm/z)'' for z < Zm and /(z) = 1 oth- 
erwise. Here, x is the fractional accuracy in the deter- 
mination of either x(-2)/'"s(zdoc) (transversal direction) or 
cff(z)-^/rs(zdcc) (radial direction) which can be obtained 
by a spectroscopic survey of volume V , measuring a galaxy 
density n at redshift z. In the above equation, D{z) is the 
growth factor, (Vb,zo,6o) are the values for a reference sur- 
vey while {xo 7) are fitted param eters obtained via 
a simulation study bv lBlake et al.l l|2006l ). which depend on 
whether one is considering a meeisurement of the acoustic 
scale in the radial or t angential direction. We employ the val- 
ues given in Table 1 o f iBlake et all (|2006l ) for a spectroscopic 
survey, as appropriate for WFMOS and the SKA. In Eq. (fTTjl 
we recognise a term oc l/\/V representing the scaling of 
the number of available Fourier modes with volume, a term 
oc l/(nD^) representing shot noise and a redshift -dependent 
cut-off term oc 1/z^ below Zm = 1.4 that suppresses non- 
linear modes (which however might also be included in 
a full non-linear analysis, thereby considerably increasing 
the BAO constraining power, see e.g. ICrocce fc Scoccimarrol 
(2008)). 

The WFMOS parameters are t aken from the results o f 
the detail ed optimization study by [Parkinson et al.l (|2007h . 
Although lParkinson et al.l (|2007l ) optimized WMFOS exper- 
imental parameters for dark energy constraints in a flat 
Universe, we expect that their general preference for a low 
redshift bin with as large as possible an area would still 
hold true even in an optimization scenario where curva- 
ture is allowed to vary. For dcfinit eness, we adopt t he va l- 
ues given m Table 2, column B of lParkinson et al.l (|20o3). 
This gives a wide bin at low redshift, covering an area of 
^low ~ 5600 deg^ at a median redshift ziow ~ 1.08, a red- 
shift width Aziow = 0.35 and a number density of galaxies 
niow = 7.1 X 10'' (/i^/Mpc^). The high-redshift bin has pa- 
rameters ^high ~ 150 deg^, Zhigh ~ 3.15, Azhigh = 0.13 
and niow = 0.13 x lO** (/i^/Mpc^). We have found that es- 
sentially all of the constraining power of this configuration 
comes from the z ~ 1 bin, in agreement with the results of 
other studies. 

The SKA is still in the design phase, hence its pre- 
cise performBMeis_somewhat uncertain at the moment (see 
e.g. lBlake et al.l (|2004l ) for an overview). We choose to repre- 
sent its capabilities by assuming measurements of both the 
transverse and radial acoustic scales equally spaced in 4 red- 
shift bins at z = 1, 2, 3, 4, each of width Az = 0.4. We further 
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assume that the SKA wiU survey the whole sky {A = 20, 000 
deg^) and that the density of galaxies will be large enough as 
to be able to neglect the shot noise term (i.e., nP > 3, where 
P is the power of the fluctuations). Some of these choices 
are somewhat optimistic and further detailed modeling is re- 
quired in order to be able to verify the capability of the SKA 
to achieve these specifications. However, we have taken here 
the SKA to represent a sort of 'ultimate' BAO measurement, 
which provides with a flavour of what the ultimate level of 
accuracy of the method might be. We add everywhere an 
external Gaussian constraint on the value of the Hubble pa- 
rameter corresponding to the HST Key Project det ermina- 
tion, i/o = 72 ± 8 km/s/Mpc (jPreedman et aLlbOOll ). 

Of course we are only dealing with statistical uncer- 
tainties here, and the issue of systematics will at some point 
have to be addressed in detail, as the statistical error be- 
comes smaller. However, BAO are particularly promising in 
this respect, thanks to the ve r y low level of systematics ex- 
pected re.g.. lTrotta fc Boweij (2006h l. 



InBoi 



Odds 



Strength of evidence 



3 CURVATURE AND BAYESIAN MODEL 
COMPARISON 

Determining whether the Universe is flat or not is one of 
the most interesting questions in modern cosmology. This is 
however not a problem of parameter constraints, but rather 
of model comparison. In this section, we briefly describe 
Bayesian model comparison and its use in forecast ing the 

J iower of future observation s (for more details, see e.g. iTrottal 
20081 ): iTrotta et all l|2008l )). We then discuss the choice of 
priors on 57^ and motivate it in the light of theoretical con- 
siderations. 



3.1 Model comparison 

From Bayes' theorem, the probability of model M given the 
data, p(A4|d), is related to the Bayesian evidence (or model 
likelihood) p{d\M) by 



p{M\d) 



p{d\M)p{M) 
p(d) 



(12) 



where p{A4) is the prior belief in model M, p{d) = 
pid\Mi)p{A4i) is a normalization constant and 



p{d\M^)^ dep{d\e,Mi)p{M^ 



(13) 



is the Bayesian evidence. Given two competing models 
M\,M2, the Bayes factor _Boi is the ratio of the models' 
evidence 



Bo 



p{d\Mo) 
p{d\Mi) 



(14) 



where large values of Boi denote a preference for Mo, 
whereas small values of _Boi denote a preference for Mi. The 
'Jeffreys' scale' (Table [5J gives an empirical scale for trans- 
lating the values of In B o\ into strengths o f belief (following 
the prescription given in iGordon fc Trottal (|2007| ) for denot- 
ing the different levels of evidence) . Recently, the framework 
of model comparison has been extend ed to include the poss i- 
bility of 'unknown models' discovery l|Starkman et al.ll2008h . 
Given two or more models, it is straightforward 



< 1.0 
1.0 
2.5 
5.0 



< 3 : 1 
~ 3 : 1 
~ 12 : 1 
~ 150 : 1 



Inconclusive 
Weak evidence 
Moderate evidence 
Strong evidence 



Table 2. Empirical scale for evaluating the strength of evidence 
when comparing two models. Mo versus JV[\ (so-called 'Jeffreys' 
scale'). The right-most column gives our convention for denoting 
the different levels of evidence above these thresholds. 



(although often computationally challenging) to com- 
pute the Bayes factor. Several numerical algorithms 
are available today to compute the Bayesian evi- 
dence. Rec ently, a v ery e ff ective algorithm, called 'nested 
sampling' (ISkflhnd |2004| . l2006l l. has become available, 
which ha s been implemented i n the cosmologi ca l con - 
text by [ Bassctt et al.' (2004*); Mukheriec ct al. (200^ 
IShaw et al., (2007); Feroz fc Hobson t2008); B ridges ot ahl 
(|2007^ . Here we are interested in the simpler scenario where 
the two models are nested, i.e. where the more complicated 
model reduces to the simpler one for a specific choice of 
the extra parameter. In our case, the extra parameter is the 
curvature, with a curved Universe reverting to a flat 
one for Jl^ = 0. Writing for the extended model parameters 
9 = (?/;, Q,K.), where the simpler (flat) model Mo is obtained 
by setting Q.^ ~ 0, and assuming further that the prior is 
separable (which is the case here), i.e. that 



p{i^,n^\Mi) = p{n^\Mi)p{i^\Mo), 

the Bayes factor can be written in all generality as 
p{Q^\d,Mi) 



Boi 



pin^lMi 



(15) 



(16) 



n„=o 



This express ion is known as the Savage-Dick ey density ratio 
(SDDR, see IVerdinefli fc WassermanI (|l995|l_ and references 
therein. For cosmological applications, see rrrottal (|2007al ')'). 
The numerator is simply the marginal posterior for f2„, eval- 
uated at the flat Universe value, Qk = 0, while the denomi- 
nator is the prior density for the model with 7^ 0, evalu- 
ated at the same point. This technique is particularly useful 
when testing for one extra parameter at a time, because 
then the marginal posterior p{QK,\d,Mi) is a 1-dimensional 
function and normalizing it to unity probability content only 
requires a 1-dimensional integral, which is computationally 
simple to do. 



3.2 A three— way model comparison 

We consider each possible choice of the curvature parameter 
K as defining a separate model. This means that we perform 
a three-way model comparison between a flat {n — 0), an 
open (k = —1) and a closed (k = -1-1) Universe. It is obvious 
that we might want to distinguish between a flat Universe 
and non-flat alternatives. However, it is also convenient to 
separate the positive and negative curvature scenarios as 
two different models. This will allow us to make statements 
on the probability that the Universe is flnite (corresponding 
to the closed case), and also to consider in a natural way a 
prior on that is flat in the log of the curvature parameter, 
as motivated below. 

For the prior probability assigned to each of the three 
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possible geometries, we make a non-committal choice of 
assigning equal probabilities to each, i.e. p{Mi) = 1/3 
(i — —1, 0, +1), where the labels of the models gives in each 
case the value of k. Of course, different choices are possible: 
for example, if one feels that inflation strongly motivates an 
almost flat Universe, this might be reflected by increasing 
the value of p{A4o) (we comment further on this below). It 
is straightforward to include such a theoretical preference 
by recalibrating our results if one wanted to. 

From the definition of the model's posterior probability, 
Eq. (|12|l and as a consequence of our assumption of equal 
prior probabilities for our models, we obtain for the posterior 
probability of the flat model the handy expression 

i + -OqI + ^0-1 

The posterior probabilities of the k 7^ models can easily 
be obtained by suitably exchanging the indexes of the Bayes 
factors. 

Each one of the models is described by a 6 parameter 
vanilla ACDM model (or a 7-parameter dark energy model 
with WeS 7^ — !)• In principle we need to specify the pri- 
ors on these parameters, too, but since they are common 
parameters to all models, their priors efl'ectively cancel, as 
shown above by Eq. (|16p . Whenever we include the extra 
parameter 7^ —1, we always add it to all models at the 
same time, therefore the model comparison is always only 
about the curvature. 

Model selection relies on a choice of prior for the extra 
parameter in the more complex model, which controls the 
strength of the Occam's razor effect, in our case f2„. Such a 
choice should be motivated by physical considerations, ide- 
ally stemming from the theoret i cal pr operties of the model 
under scrutiny (see lEfstathioul (2008) for a critical view). 
We therefore need to consider carefully the prior distribu- 
tion for the value of the parameter describing the curvature 
of spatial sections for the non-flat models. 

3.3 Priors on the curvature parameter 

A possible parameterization of the spatial curvature is given 
by the curvature parameter today, Eq. ([4]). A flat Universe 
{Hk. = 0) would therefore appear to be a point null hypoth- 
esis, to be tested against a more complex alternative model 
(with 7^ 0). In the context of inflation, however, the ge- 
ometry needs not be exactly flat. Indeed, the whole point of 
inflation is precisely to provide a mechanism to avoid such 
an implausible fine tuning. For k 7^ 0, infiation ensures that 
the curvature scale tends to zero: 

^ cxp(-2iV6) (18) 

where Nb is the number of e-folds before our c urrent co- 
movin g Hubble volume exited the horizon (see e.g. IXegmar^ 
1)2005! ) \ If we had a measure for the parameter space of 
infiationary potentials (for example from string theory), we 
could in principle convert the probability distribution for the 
potential into a prior on Nb, and from here into a prior on 
Qk- This is not necessary in practice, because local fluctu- 
ations in our Hubble volume limit the precision to which 
we can observe dev i ations from fi^ = to ~ 10~^ (see 
IWaterhouse fc ZibinI ()2008h for a more rigorous motivation 
of this result). Therefore, provided Nb > 5.8, inflationary 



predictions are observationally indistinguishable from a flat 
Universe. Given that A*';, could be anything between and 
CX3, it appears to be a reas onable approxirn ation to neglect 
models with Nb < 5.8 (see iTegmarkI ||2005l ) for a justifica- 
tion) , although such cases could be considered as a particular 
class of models if one wanted. For definiteness, in the follow- 
ing we will take the infiationary prediction to correspond to 
SIkI < 10^^, thereby extending the point null hypothesis 
that K = to include such small values of the curvature 
parameter. Because of the fundamental limitation of cosmic 
variance, we argue that it is pointless to consider the prior 
distribution of below the threshold value of 10~^. 

In summary, we describe a generic infiationary predic- 
tion as being I^kI < 10~^ (with no free parameters) and 
a prior model probability p{A4o) = 1/3. The latter assign- 
ment could of course be amended if one felt that infiation 
is compellingly motivated by its ability to solve other prob- 
lems such as the homogeneity and monopole problems, in 
which case the prior probabilities for non-inflationary mod- 
els would have to be correspondingly reducecif|. However, 
this is not essential for what follows, as we will mostly quote 
Bayes factors which give the change in degree of relative be- 
lief between two models in the light of the data. This means 
that the model's prior specification has no influence on the 
Bayes factor. In any case, it is straightforward to propagate 
a change in the models' prior probability to the model pos- 
terior probabilities that we give below. 

The model comparison is then fully defined once we 
choose a prior pdf for the extra parameter in the curved 
models, for values > 10~^. The prior should reflect our 
state of belief on the possible values of the relevant param- 
eter before we see the data. We adopt two different priors 
choices for deviations from fiatness, representing two differ- 
ent states of beliefs about the locus of possibilities for the 
geometry of the Universe. 

• Flat prior on Qk.: the 'Astronomer's prior'. This 
is prior is motivated by considerations of consistency with 
mildly informative observations on the properties of the 
Universe. Requiring that the Universe is not empty gives 
SIk > —1, barring the exotic case of a negative cosmologi- 
cal constant. The age of a Universe containing only matter 
can be approximated hytoHo = {l+fl^;^/2)~^ , which means 
that a positively curved Universe is increasingly at odds with 
the age of the oldest objects, requiring to > 10 Gyr. A pos- 
itive cosmological constant helps to solve the age problem, 
but if f2tot > 2 then to < 8/i Gyr even in a de Sitter Uni- 
verse. So unless /i ^ 1 a Universe with n„ > 1 is too young 
even in the presence of A. The lower limit for the curvature 
parameter is given by I^kI = 10""" as discussed above. How- 
ever on a linear scale this is effectively equivalent to setting 
the lower limit to 0. These considerations therefore lead to 



^ An important point is that we are here neglecting the pos- 
sibility of inflationary models predicting, for example, closed 
Uni verses with sizeab l e valu es of the curvature parameter (see 
e.g. iLasenbv fc DoranI 1I2OO5I ') for such a model). So what we de- 
scribe as a generic inflationary prediction is really only a subclass 
of possible inflationary scenarios. It would be simple to extend 
the model comparison to include other subclasses of inflationary 
models if one wanted to. 
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the prior choice: 



PAifl^Mi) = 1 for 
PA{n^\M-i) = 1 for 



-1 < n„ ^ 

^ s; 1, 



(19) 
(20) 



where the subscript A denotes that this prior is based on 
the astronomical considerations sketched above. 

• Flat prior on Infis: the 'Curvature scale prior'. 

Alternatively, we might consider the curvature scale today: 



ao 



c 



1/2 



(21) 



Clearly, a flat prior on Q.^, does not correspond to a flat prior 
on ao, as the two pdf's are related by 



p(ao) 



(22) 



Hence a flat prior on 0.^ gives an informative prior on the 
curvature scale, p(ao) oc flg''^, which prefers more strongly 
curved Universes. A flat prior on In n„ represents a state 
of belief which is indifferent with respect to the order of 
magnitude of the curvature parameter. It is easy to see that 
this implies a similar state of indifference on the order of 
magnitude of the curvature scale, since a flat prior on In 0,^ 
is flat on Inao, as well. Furthermore, such a prior is also flat 
in the number of e-folds, as a consequence of Eq. p8|l . The 
upper cutoff for the prior can be established by requiring 
that the curvature scale be larger than the Hubble horizon 
radius, H^^ . Furthermore, we are free to choose the basis in 
which the logarithm is taken, and in the following we shall 
employ base 10 logarithms. We thus define the variable 

Ok = log^Q jri^l. (23) 

These considerations lead to the prior choice 

Pc{on\Mn) = 1/5 for - 5 s; o« < (24) 

for AC = —1, 1 and where the subscript C denotes that this 
prior is based on a state of indifference with respect to the 
curvature scale. 

When employing the SDDR to evaluate the Bayes factor 
between a flat and a curved model for the curvature scale 
prior, we evaluate the marginal posterior and the prior of 
Eq. (|16|l at the value = —5 (corresponding to = 
10~^), since this is the value at which the curved models 
revert to a flat Universe for our choice of priors. 

Other choices of parameterization for curvature (and 
the associated priors) are certainly possible and might be 
well motivated from a theoretical point of view. For exam- 
ple, Adlcr & Ovorduin (2005^ introduce a constant flatness 
parameter e given by the ratio of two fundamental constants 
determining the dynamics of the expansion, and show that 
the value of e is in many ways a better indicator of "fine tun- 
ing" than Again, if one had access to the distributional 
properties of the fundamental constants and from there to 
the distribution of e, one could imagine building a phys- 
ically motivated prior on that quantity instead. However, 
since presently we are unable to predict from first principles 
the distributional properties of such quantities, we prefer to 
adopt a semi-phenomenological approach, informed by the 
physical reasoning sketched above. 



3.4 Implications for the number of Hubble 
spheres and the size of the Universe 

For closed Universes (i.e., for Mi), it is interesting to 
translate the probability distribution for Q.^ or Ok into the 
corresponding posterior for the number of particle hori- 
zon volumes_that fit int o the current spatial slice. Follow- 
ing [ScotlA^Zibi^ l|2006l '). we thus define 



Nu 



2-K 



2x - sin(2x) 



(25) 



as the ratio of the present volume of the spatial slice to the 
apparent particle horizon (assuming radiation domination 
into the infinite past), where x is the comoving radial dis- 
tance defined in Eq. ((Sjl, and for closed models G ^ x ^ ''^ ■ 
Given our choice of priors, we can easily translate the results 
of the previous section into the posterior for Nu- Clearly, un- 
der either the flat or open models {Mo ox M-i), the volume 
of the spatial slice is inflnite and hence Nu goes to infinity. 
In the Bayesian framework, we can give the probability that 
this is the case, namely that we live in an infinite Universe. 
For our choice of model priors, it follows that 

p{Nu = oo\d)=p{Mo\d)+p{M-i\d) (26) 

1 



■p{MQ\d) 1 + 



(27) 



where p{Mo\d) is given by Eq. (I17|l . For other choices of 
model priors (for example, p[Mo) ^ 1/3, representing a 
stronger degree of theoretical prejudice in favour of inflation) 
one should rescale the posteriors accordingly. 



3.5 Model comparison forecasting 

When considering the capability of future probes it is cus- 
tomary to quantify their expected performance in terms of 
a 'Figure of merit' (FOM). Several FOMs exist, but they 
mostly focus on the para meter co n strain t capabilities of fu- 
ture observations (e.g., iBassettI l|2005l )V However, many 
(and indeed perhaps most) questions of interest are actually 
about model comparison: for example, determining whether 
dark energy is a cosmological constant, or whether the Uni- 
verse is flat are clearly model comparison problems. FOMs 
geared for parameter constraints capabilities do not nec- 
essarily re flect the model com parison potential of a future 
probe (see lLiddle et al.1 (j2007l ) for details). 

A few techniques have been put forward to ass ess the 
model comparison capability of future observations: iTrottal 
(j2007bl ) has introduced a technique called PPOD, which 
computes the probability distribution o f the outcom e of a 
future model comparison; iPahud et al.l (|2006l . |2007| ) have 
looked at the ability of Planck to obtain a d ecisive model 
selecti on result regarding the spectral index; iLiddle et al.l 
|20^) have applied a similar technique to the problem of 
distinguishing between an evolving dark energy and a cos- 
mological constant. 

Here we a dopt a procedure similar in spirit 
to iLiddle et al.l (|2006t ). We want to quantify the abil- 
ity of future CMB and BAO measurements to obtain a 
correct model selection outcome about the geometry of the 
Universe. We therefore simulate data as explained above 
for three different fiducial values of Q,k'- for a flat model. 



QZ = 0, and for two different closed models, Q.* 



-10" 
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and = —5 x 10~^. From the posterior distribution 
obtained from simulated data, one can compute the cor- 
responding Bayes factor via the SDDR, Eq. (|16p . Once 
interpreted against the Jeffreys' scale, the future Bayes 
factor then allows us to determine whether the experiment 
will be accurate enough to correctly identify the true model, 
and if so with what strength of evide nce. Our pro c edure 
is th us similar to the one adopted in IPahud et al.l l|2007l . 
l2006h . 

In principle one could repeat the forecast for several 
other values of thus more densely covering the range 
of possible fiducial values. Ifowever, we found that these 3 
cases are representative of 3 interesting possibilities. The 
case Qk = —5 x 10~^ has been chosen because it lies just 
below current combined limits from CMB, BAO and SNIa, 
and within reach of the next generation of CMB and BAO 
probes. The case £7^ ~ — 10~^ is a factor of 5 below, and still 
a factor of 100 above the absolute lower limit of Q^. ~ 10^. 
Yet we will demonstrate that this scenario already presents 
very considerable challenges in terms of model confusion. 
Finally, the flat case allows us to investigate whether fu- 
ture probes can correctly determine (in a model selection 
sense) if the inflationary prediction is correct. In the follow- 
ing, we focus on the closed Universe case, because this has 
the added interest of a finite Universe, and therefore it al- 
lows to investigate the question of whether the Universe's 
spatial extent is infinite or not. In terms of parameters con- 
straints and model selection outcomes, the conclusions are 
expected to hold almost unchanged for the case of fiducial 
Universes with > 0, i.e. for the open case. 



4 RESULTS 

4.1 Current evidence for flatness 

In this section we present our model comparison analysis 
from present-day data. Our resu lts (obtained us i ng a m odi- 
fied version of the CosmoMC code. iLewis fc Bridld (|2002l l) are 
presented in Table |31 

Starting with the Astronomer's prior case, we find 
moderate evidence for a fiat Universe when compared 
with a closed model (ln_Boi ~ 4 for all cases but the 
WMAP5-I-BAO data combination with w / —1, which 
is discussed below). This corresponds to posterior odds of 
about 54:1. The evidence in favour of a fiat model is stronger 
when it is compared with the open case, as a consequence of 
the fact that the posterior for is slightly skewed towards 
values < 0, giving odds of order 200:1 in favour of the 
fiat vs the open model. When compared against each other, 
the closed model is preferred over the open model with odds 
of about 4:1. Although the odds in favour of a fiat Universe 
vs a closed one are of the same order as found in previ- 
ous works fe.g.. iKunz et al] (|2006l ) found odds of 48:1 from 
WMAP 3-yrs data and other constraints), one has to bear in 
mind that we are performing a three-way model comparison, 
while previous analyses have compared the fiat model with 
arbitrarily cur ved ones (both ope n and closed) . If we use the 
same priors as lKunz et al.l ()2006l ). we find a Bayes factor be- 
tween the fiat and curved models ln_B = 4.4 (InB — 4.6), 
for w = —1 {wcs 7^ — !)• This translates in odds of ap- 
proximately 90:1 in favour of fiatness when compared with 



a generic curved model. Thus the latest data have improved 
the model comparison outcome roughly by a factor of 2. The 
posterior probability for an infiationary, infinite Universe is 
about 98%, up from the initial 33% from our prior choice. 
The above results hold true even if one relaxes the assump- 
tion of a cosmological constant for the most constraining 
data combination, namely the one including SNIa. However, 
the evidence is favour of fiatness weakens considerably if one 
only employs WMAP5, SNIa and BAO while at the same 
time allowing for a non-constant dark energy equation of 
state (InBoi = 1.0). This is because the inclusion of BAO 
data skews the posterior for fl^ to considerably negative val- 
ues, thus preferring a closed Universe. Notice that for this 
prior the probability of a fiat Universe {p{A4o\d)) and of an 
infinite Universe ( p{Nu = oojd) ) essentially coincide, for 
the Occam's razor effect acts strongly against open models, 
as we have seen, and therefore most of the models' posterior 
probability is concentrated in the fiat Universe. 

If instead we consider the case of the Curvature scale 
prior (fiat in o„), then the Occam's razor effect penalizing 
non-fiat models is much weaker. This comes about because 
the posterior becomes fiat for Ok < —2 and stays fiat all the 
way to Ok = —5, since for such small values of the curvature 
parameter, present-day data do not provide any constraint. 
Therefore this prior choice can be seen as more conserva- 
tive in that it presents a reduced Occam's razor penalty for 
non-fiat models. From the results in Table [3] we see that for 
this prior choice the preference for fiatness is much reduced, 
although ln_Boi remains mostly positive, thus signaling a 
preference for the fiat case. For example, the odds in favour 
of fiatness vs closed (open) models are reduced to the or- 
der 3:2 (9:5), barring the case of u) ^ -1 and WMAP5, 
SNIa and BAO. However, these values are now below even 
the 'weak evidence' threshold, and therefore the model com- 
parison is inconclusive with this prior. Indeed, the posterior 
probability for the infiationary model (i.e., A4o) is now only 
about 45% (up from 33% from the prior), while the prob- 
ability of us living in an infinite Universe remains almost 
unchanged at 69% (from about 67% in the prior). This hap- 
pens because in the light of the data, the models' probability 
is redistributed in such a way that the sum of the fiat and 
open models' probability remains almost constant, despite 
the fact that the fiat model's probability has risen and the 
open model's probability has been reduced (down to about 
24% from the initial 33%). 

4.2 Constraints on the number of Hubble spheres 

For values Nu < oo (i.e., for closed models), the posterior 
probability distribution p{Nu\d, A4i) is shown in Fig. [T]for 
both choices of priors as a function of \ogj^Q{Nu)- Within the 
class of closed models, we read off Fig. [l] that the number 
of Hubble spheres is constrained to lie below Nu ~ 10® 
for the Astronomer's prior and less than Nu ~ 10^ for the 
Curvature scale prior. The sharp drop in the probability 
density for large values of Nu is a refiection of the lower 
cut-off value chosen for the priors, \ flK.\ > 10~^, while the 
difference in the upper limit is a consequence of the different 
volume of parameter space enclosed by the two priors. The 
exact value of the 99% lower limit slightly depends on the 
prior, as different priors allocate a different probability mass 
to low curvature, i.e. to large Nu - For the Astronomer's prior 
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Data sets and models 


InSoi 


InSo-i 


p{Mo\d) 


p(Nu = oo|d) 


Notes 




Astronomer's prior (flat in Q^) 


WMAP5+BAO (to = -1) 






4.1 


5.3 


0.98 


0.98 


Moderate evidence for a flat, infinite Universe 


WMAP5+BAO+SNIa (w 




-1) 


4.2 


5.3 


0.98 


0.98 


Moderate evidence for a flat, infinite Universe 


WMAP5+BAO (to ^ -1) 






1.0 


6.1 


0.74 


0.74 


Weak evidence for flatness 


WMAP5+BAO+SNIa {w 


7^ 


-1) 


3.9 


5.3 


0.98 


0.98 


Moderate evidence for flatness 




Curvature scale prior (flat in Ok) 


WMAP5+BAO (to = -1) 






0.4 


0.6 


0.45 


0.69 


Inconclusive 


WMAP5+BAO+SNIa (w 




-1) 


0.4 


0.6 


0.45 


0.69 


Inconclusive 


WMAP5+BAO (to ^ -1) 






-0.8 


0.5 


0.26 


0.42 


Inconclusive 


WMAP5+BAO+SNIa {w 


7^ 


-1) 


0.3 


0.6 


0.44 


0.67 


Inconclusive 



Table 3. Outcome of a three— way Bayesian model selection for the curvature of the Universe from current data and two choices of 
priors. For a prior choice motivated by astronomical considerations ('Astronomer's prior'), the posterior probability for a flat, infinite 
inflationary model {p{Mo\d) column) increases from the initial 33% to about 98% for the most constraining data combination, even if the 
assumption of a cosmological constant is dropped. On the contrary, the 'Curvature scale prior' returns an inconclusive model comparison, 
because in this case the Occam's razor effect is much reduced. The column p{Njj = oo\d) gives the probability of the Universe being 
infinite. 
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Figure 1. Posterior probability distribution (normalized to the 
peak) for the number of Hubble spheres contained in a spatial 
slice (for a closed Universe) from present— day CMB+BAO+SNIa 
data, for the more conservative case tOefj 7^—1 and assuming the 
Astronomer's prior (solid) or the Curvature scale prior (dashed). 



we find a 99% lower limit (1-tail) Nu > 4.8, while for tiie 
Curvature scale prior this slightly increases to Nu > 6.2 
(botli figures for the more conservative case where Wctt = — 1. 
So we conclude that, at the 99% level, the value Nu > 5 can 
be taken to be a robust lower limit to tiie number of Hubble 
spheres in the Universe. This is in good agree ment with the 
result s of the simpler analysis presented in IScott &: ZibinI 
l|2006ll , which estimated Nu > 10. 

One could also report model-averaged constraints on 
Nu, by taking into account the spread of posterior proba- 
bility between the three models: 



p{Nu\d) = pliMi\d)p{Nu\d,Mi 



(28) 



where the closed model probability, p{A4i\d), can be com- 
puted from the Bayes factors reported in Table [S] using the 



PiMi\d) = 



1 + Boi + 



(29) 



For the Astronomer's prior we obtain p(A4i\d) — 0.02 while 
for the Curvature scale prior p{Mi\d) = 0.35 for the most 
constraining data combination (and allowing for w ^ — !)• 
Because the bulk of the model probability lies with the mod- 
els where the Universe is infinite, we expect that model- 
averaged lower limits on Nu would be more stringent than 
the robust limit we reported above, but also more strongly 
prior dependent. (A similar effect is observed for model- 
av eraged constr a ints o n the dark energy equation of state 
by iLiddle et aP |200i)). For this reason, we prefer not to 
report model-averaged limits in this case. 



5 FUTURE PROSPECTS 

We now turn to the investigation of the accuracy that fu- 
ture CMB and BAO probes will achieve on Q^, both from 
the point of view of parameter constraints and, crucially, 
from the model selection perspective. Many studies have re- 
cently evaluated observational prospects using a var iety of 
probes (|Knox et al.ll2006l : lKnoxl2006l : lMao et al.ll2008l ). Here 
we improve on past works by analysing the results from a 
Bayesian model comparison viewpoint. 

We assume 3 diflferent fiducial values for Qk. '■ a fiat Uni- 
verse = 0) and two possibilities for a closed Universe, 
namely f2* = —10""^ (about one order of magnitude below 
current constraints) and a more optimistic Q% — —5 x 10"''. 
We then simulate future CMB and BAO observations as de- 
scribed above. An important point is that we simulate data 
around the true value of the fiducial model's parameters. 
This is consistent with what one would expect to obtain 
from the average of many data realizations, and analogous 
to what is usually assumed with Fisher Matrix forecasts. 
However, from the point of view of performance prediction 
and model comparison it is important to stress that this 
choice is optimistic, in that it ignores the extra uncertainty 
due to the realization noise of the specific data realization 
that one happens to observe. 
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5.1 How flat can you get? 

We first focus on tiie flat fiducial model, in which case no 
deviation from flatness should be observed (with the impor- 
tant caveat of realization uncertainty given above) and fu- 
ture probes will further tighten constraints around Q^. = 0. 
In Table H] we report the projected posterior la constraint 
on r2„ as well as the 99% (2.58ct) 1-tail lower limit on 
This quantity would be the appropriate figure to report in 
the case that no deviation from flatness is found and one 
wanted to constrain positively closed models at the 99% 
level. This limit can also be translated into the correspond- 
ing 99% lower bound on the number of Hubble spheres, Nu, 
which is also given in Table U Combination of future CMB 
data with WFMOS BAO determinations will constrain cur- 
vature at the ~ IQ-^ level, with the degradation in accuracy 
coming from dropping the assumption of a cosmological con- 
stant being about a factor of 2. Interestingly, once Planck 
data are available there is not much to be gained in terms 
of curvature constraints from a CVL CMB experiment. An 
SKA-like BAO experiment will further tighten constraints 
by a factor of about 5, and reduce the dependency of the 
marginal curvature accuracy on the assumptions about the 
dark energy equation of state. Constraints in the {Qk.,Wcs) 
plane for the flat prior case are depicted in Fig. O showing 
how SKA will essentially eliminate the correlation between 
the two parameters, leading to independent constraints on 
the curvature and the effective dark energy equation of state. 

This result could potentially be weakened if one allowed 
for a more general dark energy time dependence than we 
have considered here. However, iKnox et al.l (|2006l ) showed 
that WFMOS BAO constraints on Jl^ are remarkably ro- 
bust even if one allows for an evolving dark energy of the 
form w{z) = Wo + Waj^- This is mainly due to the extra 
constraining power coming from the high-redshift bin, which 
in our analysis played a subdominant role since we assumed 
that the effective equation of state is constant with red- 
shift. Even f or the more general {wo,Wa) parameterization, 
iKnox et al] l|2006h found that WFMOS-like constraints on 
Qk. are only degraded by about 50% wrt our result (see be- 
low for further comments about the impact on model confu- 
sion). We note that our forecast for SKA-like BAO data is 
of the same order of the accuracy that could be achieved by 
a combination of weak lensing and BAO observations by the 
Large Synoptic Survey Telescope (LSST) whe n marginaliz- 
ing over a more general {wo,Wa) jZhanl 120061 ). However, if 
one models the dark energy equation of state as a contin- 
uous function, then constraints on curvature are very con- 
siderably degraded. Even a combination of weak lensing and 
BAO observations by the LSST w ill only achieve a relatively 
modest accuracy ~ 0.017 on (|Zhan et al.ll2009l '). 

In terms of constraining the number of Hubble spheres, 
WFMOS data would increase the current lower limit Nu > 5 
by almost 2 orders of magnitude to Nu > 300 — 400, while 
SKA-like BAO observations would further improve this by 
1 order of magnitude to Nu > 2000 — 3000 (all figures are 
given before model-averaging). 

We now turn to the model selection question of whether 
future experiments will be able to determine unambiguously 
that the Universe is flat, should this be the case. Results are 
shown in Table O which gives value of In Boi , the Bayes fac- 
tor between the (correct) flat model and the closed model 
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CVL + WFMOS 



o 
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o 
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Figure 2. Future constraints on curvature and the dark energy 
effective equation of state from various combinations of future 
probes, for the Astronomer's prior. Contours delimit 68% and 
95% joint credible regions, the cross gives the fiducial value. 
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Planck 


+ WFMOS 


1.76 
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-3 


-4.17 




392 


CVL -1- 


WFMOS 


1.60 


10 


-3 


-3.85 


10-3 


443 


Planck 


+ SKA 


5.64 


10 


-4 


-1.34 


10-3 


1970 


CVL -1- 


SKA 


4.58 


10 


-4 


-1.07 


10-3 


2732 






Planck 


-1- WFMOS 


2.22 


10" 


-3 


-4.58 


10-3 


284 


CVL + 


WFMOS 


2.08 


10" 


-3 


-4.40 


10-3 


293 


Planck 


+ SKA 


6.38 


10 


-4 


-1.50 


10-3 


1676 


CVL + 


SKA 


4.58 


10 


-4 


-1.05 


10-3 


2723 



Table 4. Posterior constraints on the curvature parameter 
from future CMB and BAO probes, taking a fiducial value O* = 
and for a fiat prior on 51^. 



(recall that InBoi > favours the flat model). The values 
of InBo-i are within a few percent from the ones given in 
the table, and hence are not displayed (the small difference 
comes from the fact that the dependency of the observables 
is not precisely symmetric in f2„ around 0,^ ~ 0)- Our find- 
ings show that all of the experiments will be able to return 
strong evidence (In Boi > 5) for the case of a fiat prior on 
JIk- This results holds true even if we relax the assumption 
that dark energy is in the form of a cosmological constant. 

However, the strength of evidence is much reduced if 
instead one employs a prior that is flat on o^, as shown 
in the right-hand-side of Table [S] ('Curvature scale prior'). 
Even the most constraining experiments (CVL-I-SKA) will 
struggle to gather weak evidence (InBoi > 1.0) in favour of 
flatness. This comes about for two reasons. First, evidence 
accumulates only proportionally to the inverse error on the 
parameter of interest, hence in the Bayesian framework it is 
much easier to disprove a model (where the evidence goes 
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(n* = 0.0) 

Probe 

Planck + WFMOS 
CVL + WFMOS 
Planck + SKA 
CVL + SKA 



Astronomer's prior 
ui = —1 w —1 

6^0 5^9 

6.2 5.9 

7.1 6.2 

7.5 6.3 



Curvature scale prior 
w = —1 w ^ —1 

OT 07 

0.8 0.7 

1.0 1.0 

1.1 1.1 



{n% = -5 X 10-3) 

Probe 

Planck + WFMOS 
CVL + WFMOS 
Planck + SKA 
CVL + SKA 



Astronomer's prior 
UI = —1 w 7^ —1 

T3 2^6 

0.4 2.0 

-34 -22 

-55 -50 



Curvature scale prior 
w = —1 w 7^ —1 

-3.5 ~n 

-4.5 -2.0 
-50 -40 
-65 -58 



Table 5. Outcome of Bayesian model selection from future data, 
generated from a flat Universe. The table gives values of In Bqi , 
the Bayes factor between a flat and a closed model, using a flat 
prior on Qk ('Astronomer's prior') or a flat prior on ('Curva- 
ture scale prior'). The analysis with a flat prior on Qk gives strong 
evidence in favour of the flat model even when the assumption of 
a cosmological constant is relaxed (ui 7^ — 1 column), while using 
a flat prior on Ok the strength of evidence is just above the 'weak' 
threshold even for the most powerful probe (CVL+SKA). 



o 



0.6 



0.4 - 



0.2 - 








-4 -3 
Logio(^J 



Figure 3. Normalized posterior on assuming a fiducial value 
O* = 0, reconstructed using a flat prior on 0^ and assuming w = 
— 1, for different combinations of future data. From right to left: 
Planck+ WFMOS (solid), CVL+WFMOS (dotted), Planck+SKA 
(short-dashed), CVL+SKA (long-dashed). 



exponentially in the number of sigma discrepancy with the 
prediction) than to confirm it. Secondly, the Occam's razor 
effect penaUsing non-flat models is much reduced under the 
assumption of a flat prior in Ok, as the net result of two 
opposite effects. Fig. [3] shows the posterior pdf on o„ for 
the different probes (assuming w — —1). It is clear that for 
values of <ti —3 the posterior becomes essentially flat, 
reflecting the inability of the experiment to measure a cur- 
vature value much below that threshold. At the same time, 
the volume of parameter space enclosed by a prior flat in log 
space is increased with respect to the case of a linear scale. 
This ought to favour the simpler (flat) model. But the pos- 
terior volume is also increased, and therefore the net effect 
is to reduce the overall Occam's razor penalty term (which 
goes as the log of the ratio between the two volumes), hence 
the strength of evidence in favour of flatness is reduced. 



Table 6. Outcome of Bayesian model selection from future data, 
generated from a closed Universe, f2* = —5 X 10^'^. The table 
gives values of InBoi, the Bayes factor between a flat and a 
closed Universe. Negative values correctly favour the closed case, 
while positive values wrongly favour the flat case, giving rise to 
model confusion. SKA-quality BAO data are required to over- 
come model confusion independently of the choice of prior. 

5.2 The danger of model confusion 

We now turn to the case where the fiducial model is closed, 
and evaluate the resulting evidence from future data. In this 
case, a successful model comparison should return a prefer- 
ence for the closed model. 

We start with the more optimistic case of a relatively 
large fiducial value for the curvature parameter, f2* — 
—5 X 10~^, roughly a factor of 2 below present-day con- 
straints. We give results in Table |6] which shows that the 
flat prior on Ok. always returns the correct model comparison 
(negative values of InBoi in the table). However, the As- 
tronomer's prior incorrectly penalizes curved models when 
the constraining power of the data is insufficient to over- 
turn the Occam's razor effect (positive values in the table). 
This 'model confusion' effect is worse when the equation 
of state of dark energy is allowed to change, in which case 
e.g. Planck -(-WFMOS would incorrectly gather moderate ev- 
idence in favour of flatness. With CMB and WFMOS data, 
the analysis is subject to model ambiguity, i.e. the result de- 
pends on the choice of prior. In order to recover the correct 
model selection outcome unambiguously, one needs SKA- 
quality BAO data to complement the CMB distance probes 
(negative values of InBoi for both priors). 

The danger of model confusion becomes stronger the 
smaller the flducial value one chooses for fi*. We illustrate 
this by considering our third fiducial model, namely a closed 
Universe with Q% = —10'^, which is about one order of 
magnitude below current limits but still two orders of mag- 
nitudes above the fundamental fluctuation limit. The model 
comparison outcome is given in Table [T] which shows that 
this case results in widespread model confusion for the As- 
tronomer's prior, for which the flat Universe is incorrectly 
preferred with moderate to strong evidence by all combina- 
tions of probes. For the Curvature scale prior, instead, the 
outcome is always inconclusive, even though the CVL-I-SKA 
combination does reach the correct conclusion, albeit with 
evidence which falls short even of the 'weak' threshold. 

Some comments are in order about the robustness 
of those results with respect to changes in the assumed 
dark energy model. In particular, an evolving dark energy 
comp onent could mimic to an extent the effect of curva- 
ture ijClarkson et al.l[2007l ). and this would lead to increased 
uncertainty in the curvature parameter and thus to increased 
model confusion. To estimate the impact of this effect, we 
have repeated the analysis for a subset of the cases discussed 
above, but marginalising over a 2-parameters dark energy 
equation of state of the form w{z) = wo + WajTr^- In this 
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Probe 


-10-3) 


Astronomer's prior 
UI = —1 w —1 


Curvature scale prior 
w = —1 w ^ —1 


Planck 


+ WFMOS 


5.6 


5.5 


0.6 


0.6 


CVL + 


WFMOS 


5.6 


5.2 


0.4 


0.6 


Planck 


+ SKA 


5.0 


5.2 


0.0 


0.1 


CVL + 


SKA 


4.4 


4.4 


-0.6 


-0.6 



Table 7. As in Table El but for a fiducial value C* = -10"^ 
For such a small value of the curvature, only the CVL+SKA 
data combination achieves the correct model selection (albeit with 
undecided odds) and this only when employing the curvature scale 
prior. All other cases are subject to model confusion. 



case, the values of In Boi are reduced by ~ 10% — 20% with 
respect to the case where a Wcs 7^—1 model was assumed. 
This change is not large enough to modify in a significant 
way the outcome of model selection reported in Tables [5HZ1 
Therefore we conclude that assuming a more general dark 
energy equation of state does not impact very strongly on 
our results about the danger of model confusion. 



5.3 Avoiding model confusion 

In the hght of the findings in the previous section, it is in- 
teresting to estimate the required accuracy on f2„ in order 
to ensure that future probes will not be subject to model 
confusion. For a given fiducial value of|r2„| >10~^,we wish 
to estimate the accuracy needed so that the model com- 
parison correctly favours the closed model over a flat one 
independently of the choice of prior. 

This can be achieved by using a Gaussian approxima- 
tion to the future likelihood and employing the SDDR to es- 
timate the Bayes factor between the closed and open model 
that a future experiment would obtain. We start by consid- 
ering the case of the Astronomer's prior. We assume that the 
marginal likelihood of is approximately described by a 
Gaussian with mean Q% (i.e., centered on the fiducial valu^l) 
and variance E^. Then adopting the Astronomer's prior, the 
Bayes factor between the flat and the closed model is given 
by, from Eq. ([T6)l. 



InSoi 



E 1 i^iy 



(30) 



where AO = 1 is the width of the Astronomer's prior on the 
curvature parameter and the last term of the right-hand- 
side is defined as 



Erf 



Ae- \0 . 
v/2E 



and Erf(a;) denotes the error function, 

Erf(a;) = - [ exp(r^)dr. 



(31) 



(32) 
" Jo 

The function /a accounts for the upper and lower limits in 

3 As discussed above, this neglects the realization noise and is 
therefore equivalent to an ensemble-averaged forecast, analogous 
to what is usually done for Fisher matrix forecasts. However, nu- 
merical investigations suggest that realization noise is a subdomi- 
nant source of uncertainty in this context (Andrew Liddle, private 
communication) . 



-1 
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Model 
confusion 
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Correct 
model 
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-3 -2 -1 
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Figure 4. Bayes factor from future data for the Astronomer's 
prior as a function of the true value of the curvature parameter 
fi* and the future marginal accuracy on Ok, E. The red, thick 
line separates regions of model confusion (above the line, shaded, 
ln_Boi > 0, wrongly favouring a fiat Universe) from regions of 
correct model selection (ln_Boi < 0, white, correctly returning a 
preference for a closed Universe). The contours denote increasing 
levels of evidence, with values of In Bqi as labelled. The contours 
below the red line delimit regions of weak, moderate and strong 
preference for the closed Universe (from top to bottom). 



the Astronomer's prior distribution when computing the ev- 
idence. It is easy to see that when the posterior is sharply 
localized within the prior, i.e. for f2*/E ^ 1 it follows that 
/A(r2^,E) iln27r. Thus in Eq. ((30}, the first two terms 
on the right-hand-side represent the Occam's razor effect 
(notice that — InE/AO > 0, thus favouring the flat model) 
, while the last term describes the relative quality of fit be- 
tween the closed and fiat model. We have checked the accu- 
racy of the approximation of Eq. (|30p against the full numer- 
ical results in Tables [6] and [7] adopting the error estimates 
given in Table and we have found it to be excellent, with 
an accuracy of a few percent. 

The result is plotted in Fig.|^ where the red, thick line is 
the contour level In Boi = which separate the region where 
the model comparison correctly favours a closed Universe 
(bottom right corner, in white) from the 'model confusion' 
region, where the flat Universe is incorrectly preferred due 
to the Occam's razor effect (shaded region above the red 
line). For a given value of the curvature parameter on the 
horizontal axis, the red line thus gives the required marginal 
accuracy on r2„ to avoid model confusion. We will come back 
below to discussing what this means in terms of the required 
discovery threshold. The unfilled contours below the red line 
denote values InBoi = —1.0, —2.5, —5.0 (weak, moderate 
and strong evidence for curvature, respectively, from top to 
bottom in the figure). Because the evidence against the null 
hypothesis of a fiat Universe grows exponentially in the tails 
of the distribution, those contours are relatively close to the 
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Figure 5. As in Fig. |4l but for tlie Curvature scale prior. The 
blue, thick line separates regions of model confusion (above the 
line) from regions of correct model selection (InSoi < 0). The 
shaded areas denote regions of increasing model confusion, from 
light to dark. 

In Boi = threshold. This means that a relatively modest 
increase in accuracy can lead to 'strong' evidence in favour 
of curvature. The situation is not symmetric with respect 
to the null hypothesis: the evidence increases only linearly 
with the accuracy in case of a null result, hence it takes a 
much larger accuracy to accumulate evidence in favour of 
the null. This is reflected by the larger spacing between the 
evidence contours in the model confusion region. 

Turning now to the case of the Curvature scale prior, 
the Bayes factor can be computed in an analogous fash- 
ion, by replacing the flat prior on £7^ by the prior equiv- 
alent to a flat prior on o„, namely piyiK) ~ M/Q.K (for 
— 1 ^ SIk ^ — 10~^) and M is a normalization constant. The 
Bayes factor can then be computed numerically using the 
SDDR, Eq. p6|l . The resulting outcome for model selection 
is shown in Fig. [5] where the blue, thick line again separates 
region of correct model selection from regions of model con- 
fusion, as a function of the fiducial value for the curvature 
and of the marginal accuracy. By comparing with Fig. HI 
we notice that the Curvature scale prior is less subject to 
model confusion than the Astronomer's prior, since for the 
former the strength of evidence in favour of a flat Universe 
is lower in the model confusion region and it barely reaches 
the 'moderate' evidence threshold. Furthermore, the blue 
fine is always above the red line (see Fig. [6]), which means 
that model confusion for the Curvature scale prior is avoided 
with less stringent requirements on the marginal accuracy 
E than for the Astronomer's prior. 

5.4 Limits to the knowability of the the geometry 

The comparison between the two priors is further inves- 
tigated in Fig. [S] where we plot the contours separating 




-5 -4 -3 -2 -1 

Figure 6. The red (blue) line delimits regions of model confu- 
sion (above the line) for the Astronomer's prior (Curvature scale 
prior). The light— shaded region between the red and blue lines is 
a zone of 'model ambiguity', where the model comparison results 
depends on the prior. The diagonal, dotted lines denote approx- 
imate regions of Sct and 5(t (from top to bottom) discovery for 
a given true value of the curvature parameter, f2* . A ba discov- 
ery threshold guarantees an ambiguity- and confusion-free model 
comparison outcome. The horizontal lines give the expected ac- 
curacy of future CMB and BAG probes. 



the model confusion region for both priors (red for the As- 
tronomer's prior and blue for the Curvature scale prior). The 
dark shaded region labeled 'model confusion' leads to an er- 
roneous model comparison result for both priors, while the 
light shaded region between the two lines is a zone of 'model 
ambiguity' — where the outcome of model comparison de- 
pends on the choice of prior. In such a case, better data (i.e., 
smaller E) are required in order to resolve the ambiguity. It 
is interesting to investigate the accuracy necessary to obtain 
an ambiguity- and confusion-free model selection. In Fig. |6l 
the diagonal, dashed lines represents approximate 3 and 5(t 
detection thresholds (from top to bottom). We can see that, 
except for fairly large values of Q.*^ > 0.1, a 3cr 'detection' is 
subject to both model ambiguity and model confusion. On 
the other hand, a 5(t detection leads to an unambiguous and 
correct model choice all the way down to f2J > 7 x 10~^. 

This is further substantiated by the results tabulated 
in Table |8] giving the required accuracy (both in absolute 
value and number of a discovery) to achieve moderate or 
strong evidence under both priors, for a few representa- 
tive choices of the fiducial curvature value. Under the As- 
tronomer's prior, moderate evidence in favour of curvature 
requires a ~ 4a detection, while for strong evidence a ~ 5(t 
detection is necessary. As mentioned above, the Curvature 
scale prior is less demanding due to its reduced Occam's ra- 
zor effect: moderate evidence is achieved with a ~ 3.2(t de- 
tection threshold, while strong evidence is obtained at ~ 4(j. 
This of course comes at the price of a much reduced evidence 
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in favour of a flat Universe if that is indeed the true model, 
as discussed in connection with the results presented in Ta- 
ble [5] 

In conclusion, our results imply that a 5a detection 
threshold ought to be recommended in order to obtain a se- 
cure and ambiguity-free model selection. It is perhaps amus- 
ing that a full Bayesian treatment of the problem concludes 
that the 5a detection threshold traditionally used in par- 
ticle physics (with its frequentist framework) ought to be 
employed. 

Finally, we can revise the conclusion about the funda- 
mental limit to the knowability of the geometry of the Uni- 
verse. It is usually argued that this is of order ISIkI ~ 10^^, 
because this is the typical size of curvature fluctuations 
due to primordial inhomogeneities. However, Fig. [S] shows 
that model confusion sets in for value of the curvature 
< 10-*, which means that if the true value of the 
curvature is below this threshold we will not be able to 
gather evidence for it. We conclude that the fundamental 
limit to our ability to detect the curvature of the Universe 
(if present) is of the order |f2„| ~ IQ-*, which is an order 
of magnitude greater than previous estimates. Below that 
value, the Occam's razor arguments inbuilt into Bayesian 
model selection imply that we ought to revert to preferring 
a fiat Universe. Therefore if the curvature is in the 'unde- 
cidable interval' 10"^ ^ < 10'" no amount of data will 
be able to determine that the Universe is non-flat. 



6 CONCLUSIONS 

We have subjected the geometry of the Universe to a detailed 
scrutiny from a model comparison perspective, performing 
a three-way model selection with two physically motivated 
priors. We found that present-day data provide up to mod- 
erate evidence in favour of flatness (maximum odds of 66:1) 
for a specific choice of prior (the Astronomer's prior) and 
assuming that dark energy is a cosmological constant. A 
Curvature scale prior and a relaxation of the assumption on 
the nature of dark energy weaken this result considerably, 
giving only inconclusive odds of less than 3:2 in favour of 
flatness. Correspondingly, the probability that the Universe 
is infinite lies in the range from 67% to 98%, depending on 
assumptions. If the Universe is not infinite, we have found a 
robust lower limit to the number of Hubble spheres, Nu > 5. 

We have discussed the prospects for future CMB and 
BAO probes to determine with strong evidence the geome- 
try of the Universe. CMB data coupled with WFMOS BAO 
observations will achieve an accuracy on Sl^ of the order 
~ 1 - 2 X 10"^, while SKA-like BAO data will further in- 
crease the accuracy to ~ 4 — 6 x 10-*. Allowing for the effec- 
tive equation of state of dark energy to be different from —1 
(although constant in redshift) will not significantly decrease 
the accuracy with which CMB-I-SKA data will determine 

Finally, we have shown that a model selection perspec- 
tive places much more taxing requirements on the accuracy 
of future datasets than one would naively assume. In par- 
ticular, a 5a detection threshold is recommended in order 
to avoid both model confusion and model ambiguity in the 
determination of the geometry. However, if the value of the 
curvature parameter is smaller than ~ 10"* we found that 



no amount of observations will be able to decide on the 
true geometry of the Universe. Achieving this lower limit 
will require an improvement of another factor of 20 over 
what a CVL CMB experiment with an SKA-like BAO probe 
will obtain. This might be feasible once other, orthogonal 
datasets such as weak lensing and SNIa observations are 
added to the likelihood, although it will be a formidable 
challenge to control systematics at this level of statistical 
accuracy. 
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True value 


Astronomer's prior 
moderate evidence strong evidence 


Curvature scale prior 
moderate evidence strong evidence 
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-5 X 10"^ 

-10-3 

_10-4 


1.23 X 10-3 (4.06cr) 
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3.12 X 10-"' (3.170-) 
2.79 X 10-5 (3.22(t) 


1.26 X 10-3 (3.95cr) 
2.50 X lO-"' (3.95cr) 
2.26 X 10-5 (3.990-) 



Table 8. Required accuracy on the marginal error on to avoid model confusion and to achieve different thresholds of evidence in 
favour of a closed Universe for the two priors considered in the text (absolute value and relative number of a in parenthesis). A 5a 
discovery threshold guarantees an ambiguity- and confusion-free model comparison outcome down to |f2* | = 10—*. 
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